SHEP-0518, UW/PT 05-16 



Finite- Volume Effects for Two-Hadron States in 

Moving Frames 



C.h. Kim a , C.T. Sachrajda a and Stephen R. Sharpe 

^ ■ a School of Physics and Astronomy, Univ. of Southampton, 

<3 Southampton, S017 IB J, UK. 

(N 



a.b 



OO 



b 



Department of Physics, University of Washington, Seattle, WA-98195-1550, USA . 



> 

^ ■ Abstract 
O 

t^. ' We determine the finite-volume corrections to the spectrum and matrix elements of two- 
hadron states in a moving frame, i.e. one in which the total momentum of the two- 
q . hadrons is non-zero. The analysis is performed entirely within field theory and the results 
are accurate up to exponential corrections in the volume. Our results for the spectrum 
are equivalent to those of Rummukainen and Gottlieb which had been obtained using a 



o ,| relativistic quantum mechanical approach. A technical step in our analysis is a simple 
derivation of the summation formulae relating the loop summations over the momenta of 
the two hadrons in finite volume to the corresponding integrals in infinite volume. 

X 
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1 Introduction 



Lattice simulations are necessarily performed in a finite volume. In the evaluation of many 
quantities of phenomeno logical interest, such as hadronic masses or matrix elements with 
at most a single hadron in the external states, the errors caused by the finiteness of the 
volume decrease exponentially as the volume is increased [1]. In practice, as long as 
the size of the lattice is large enough compared to the range of the strong interactions, 
these errors are numerically negligible. This is not the case, however, when two (or 
more) hadrons are present in external states. For such states the finite-volume effects 
decrease more slowly, as powers of the box size, L, and need to be understood in order to 
obtain physical quantities with good precision. The theory of such effects has been fully 
developed for two particles in their rest frame, i.e. with total momentum P = 0. The 
spectrum of such states was worked out in refs. [2-4] and the finite volume corrections to 
the matrix elements were obtained in refs. [5,6]. 

This paper concerns the extension of these results to two-particle states with total mo- 

— * 

mentum P ^ 0, which we call a moving frame. A generalisation of the results for the 
spectrum of such states to a moving frame was proposed some time ago in ref. [7]. This 
was based upon a relativistic quantum mechanical approach in which the two particles 
satisfied a Klein-Gordon equation with an interaction potential of finite range. Part of 
the motivation for our project was an attempt to understand some of the steps in the 
derivation of ref. [7]. We have therefore revisited the problem using an entirely field the- 
oretic approach. Although our final result looks superficially quite different from that 
obtained in ref. [7], we find that it is possible to transform our expressions into theirs. 
Thus we confirm, albeit indirectly, their assumptions and methodology. Our derivation 
and results are useful for a number of reasons. First, the extension of the results for finite 
volume corrections to matrix elements are new. Second, we give an alternative and simple 
derivation of the summation formulae that are needed. And, finally, our derivation of the 
quantization condition for two-particle states is simple and transparent, and in particular 
does not require the use of an intermediate quantum mechanical theory. 

There are several reasons why it is important to have an extension of the theory of fi- 
nite volume corrections to a moving frame. Some were discussed in ref. [7]: one obtains 
additional results for the scattering phase shifts at different values of the centre-of-mass 
momenta from the same lattice configurations; one can more easily study p-wave decays 
such as p — > 7T7r; and one avoids the need for vacuum subtractions in s- waves. An ad- 
ditional important motivation is provided by the need to control non-perturbative QCD 
effects in K — > nn decays, and hence, for example, to understand quantitatively phenom- 
ena such as the AI = 1/2 rule and the value of e'/e. Whilst it is true that the decay 
amplitudes can in principle be evaluated in the rest-frame, the ability to perform simu- 
lations in a moving frame will be very important in controlling systematic uncertainties. 
For example chiral perturbation theory is an important guide in performing the extrap- 
olation of lattice results to those corresponding to physical u and d quark masses. In 
this approach there are a number of unknown constants, and results from moving frames 
provide considerably more data points to allow the determination of these low-energy con- 
stants (see, for example, ref. [8]). Furthermore, as noted above in the case of scattering 
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amplitudes, the use of a moving frame obviates the need for vacuum subtractions in the 
isoscalar channel. This is important since such subtractions are under poor statistical 
control. 

Throughout this paper we assume that the fields satisfy periodic boundary conditions in 
the moving frame (and we take the temporal extent of the lattice to be infinite). There has 
recently been discussion of the possibility of using anti-periodic [9] or twisted boundary 
conditions (if)(xi + L) = exp(i9i)ip(xi) i = 1,2,3) in lattice simulations [10]. Although 
we do not consider such boundary conditions explicitly in this paper, the techniques 
presented below, and in particular the summation formulae, can readily be generalised 
for the different cases. 

We now briefly present our result for the finite-volume corrections to the two-particle 
spectrum in a moving frame. The example we have in mind is of two pions, and so, for 
concreteness, we will call the particles "pions" throughout this paper. We note, however, 
that the results in this paper can readily be generalised to the case of two particles with 
different masses (mi and m 2 say), as long as they are below the inelastic threshold. 
An important example is the nucleon-pion system. The quantization condition given 
below (see eq. (|50|) for the general case) in terms of the centre-of-mass momentum of each 
particle, q*, remains valid, q* is then related to the total centre-of-mass energy E* by 
the standard kinematical relation 4q* 2 = E* 2 — 2{m\ + mj) + {m\ — m\) 2 / E* 2 . For 
the remainder of this paper we limit the discussion to the case when the two hadrons 
are degenerate (mi = m-i = m) but it should be remembered that the generalisation to 
non-degenerate particles is straightforward. 

The finite-volume cubic box breaks the Euclidean symmetry and hence mixes different 
partial waves. For P / 0, lattice states which are predominantly s-wave also contain 
components with I — 2, 4, . . . (in contrast to the case with P = where the mixing begins 
at I — 4). The results presented here were obtained under the assumption that only 
scattering in the s-wave is important in the energy range of interest (which is a good 
approximation for K — > txtt decays). The general results are given in the body of the 
paper. Let the two-pion state have total momentum P and let E be one of the energy 
eigenvalues in the finite- volume V = L 3 . The total energy in the centre-of-mass frame 
is denoted by E* (E* 2 = E 2 - P 2 ), and q* is defined by E* 2 = A[q* 2 + m 2 } (2q* is the 
magnitude of the relative momentum). We find that the spectrum of two-pion states is 
determined by the relation 



where S(q*) is the physical s-wave phase-shift and the function (p p (q*) is defined by 






(2) 



with the box size entering through the following regularized sum 



The summation in this expression runs over the usual finite- volume momenta for periodic 
boundary conditions: k = (2tt / L) n, with n a vector of integers. The corresponding 



pion energy is uj^ = y k 2 + m 2 . The quantities (ujl,k*) are the energy and momentum 
of the four-vector (u^., k) boosted into the centre-of-mass frame. V indicates that one 
should take the principal value of the integration and the superscript P recalls that we 
are considering a two-pion system in a moving frame. The exponential factors in eq. (|3*|) 
regulate the ultraviolet divergences in the sum and integral, and these divergences cancel 
in the difference of the two. The remaining a dependence is exponentially suppressed as 
L — > oo. Other choices of ultraviolet regulator are of course possible, e.g. the analytic 
regularization used in refs. [2,3,7]. For this presentation we prefer to keep a regulator 
such as the one in (jSJ), so that the cancellation of UV divergences between the sum and 
the integral is manifest. 

Our result for the quantization condition, eq. (0), differs from that of ref. [7] in the form 
of the kinematical function. Our c p (q* 2 ) of eq. © is replaced in ref. [7] by 

X/ "1*2 ~2 W 



7L 3 g q* 2 — r 



r2 = -lh--) + (5) 



where 7 is the Lorentz factor for the transformation between the moving and centre-of- 
mass frames, and 1 

1 /. 

with k\\ and k± the components of k parallel and perpendicular to P, respectively. The 
ultraviolet divergences in eq. (0} are regulated by analytical regularization. In sec. 0] we 
demonstrate that the expression for c p (q* 2 ) in eq. (jSJ) is equal to that in eq. (jlj), up to 
terms which vanish exponentially with the volume. The two quantization conditions are 
therefore equivalent. 

We show in sec.0]that this equivalence also extends to the general situation in which there 
are multiple partial waves with non-vanishing phase shifts. In light of this, the discussion 
given in ref. [7] concerning the group theoretic constraints on the mixing between partial 
waves, and the numerical values of the kinematical functions remains valid. Thus we do 
not repeat the discussion in this paper. 

In earlier papers tan(0 p ) was frequently written in terms of the dimensionless variable 
g*(L/27r) (and m(L/27r) where appropriate). In order to avoid confusion when comparing 
our results with earlier ones, we underline that our notation is different and the energy 
and momentum variables in this paper (such as k, u^, k*, u>l, q*, r) are all dimensionful. 2 

The other principal result of this paper is the expression for the finite- volume corrections 
in matrix elements of local composite operators between external states, one or both of 

1 Ref. [7] uses a slightly different definition of r in which the P/2 term in eq. (JSJ comes with a positive 
sign. This difference has no effect on the sum in eq. ifljl. however, since a change of the summation 
variable k —* k + P brings our form into that used in Ref. [7] . We find the form in eq. 10 more natural, 
as will become clear in the subsequent discussion. 

2 Of course this is a matter of choice. We choose to use dimensionful variables because it is the physical 
energies and momenta which we wish to keep fixed as L is increased. 
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which contain two hadrons (with an energy below the inelastic threshold). An important 
application is to K — > nn matrix elements. Our results [see eq. (|§U|) below] generalise 
those of Lellouch and Liischer [5] and Lin et al. [6] for the centre-of-mass frame. 

The issues addressed in this paper have also been considered by Christ, Kim and Ya- 
mazaki. Using a related but different approach, they have also confirmed the validity 
of the results of ref. [7], and derived the generalisation of the result of Lellouch and 
Liischer for moving frames. Their paper is being released simultaneously with the present 
article [11]. 

The remainder of this paper is organized as follows. In the next section we derive the 
summation formulae that are needed to evaluate the generic loop summation over the 
momenta of two-particle states. In sec. H3 we use the summation formulae to derive the 
quantization condition for two-pion states in finite volume. We then show, in sec. EJ 
that our results are equivalent to those obtained in ref. [7], and in sec. El we verify the 
quantization condition at lowest non-trivial order in perturbation theory in a A0 4 model 
of pion interactions. Section |U] contains our derivations of the finite volume corrections to 
K — > titi matrix elements. We present our conclusions in sec. 



2 Summation Formulae 

Key ingredients in the calculation of finite-volume effects for two-hadron states are sum- 
mation formulae relating the sums over the discrete momenta in a finite volume to the 
corresponding integrals over the continuous spectrum in infinite volume. We start from 
the Poisson summation formula, 

where the summation on the left-hand-side is over all integer values of n = (111,712,113), 
with k = (2ir/L) n, while that on the right-hand-side is over integer values of I = l 2 , 13) 
excluding I = (0,0,0). In the following, we consider functions f(k ) whose Fourier trans- 
forms, f(r), are non-singular, and are either contained in a finite spatial region or decrease 
exponentially as \r | — > 00. If we apply the Poisson summation formula to such functions, 
the terms with I 7^ on the right-hand-side of eq. © decrease at least exponentially as 
the box size is sent to infinity, so that 

up to exponentially small corrections. We note that functions f(k) with these properties 
have no singularities for real k, and fall off fast enough at \k\ ^00 that the integrals in 
eq. (JHJ) converge. 

We begin our discussion with a review of the summation formulae in the centre-of-mass 
frame and then in section 12.21 generalise the discussion to the moving frame. 
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2.1 Summation Formulae in the Centre-of-Mass Frame 



For two-hadron correlators with an energy below the inelastic threshold, the finite- volume 
corrections are contained in summations of the form 3 

k H 

where we assume that q 2 is such that there is no term in the sum with k 2 = | k \ 2 = q 2 
and that f(k ) has the properties discussed above, and in particular has no singularities 
for real k. The singularity at k 2 = q 2 forbids the simple replacement of the sum with an 
integral. Defining spherical coordinates, k = (k, 8, <p), we expand / in spherical harmonics 

oo / 

= £ E fUk)k l Y lm (e,4>), (9) 

1=0 m=-l 

so that (with q = \q\) 

S(q) = E Si m (q) where S lm (q) = 4 E k * Y iM <P) ■ W 

l,m k 1 

In the following, it is convenient to consider each Si m (q) individually. 

To make use of eq. (J7J), we subtract from the summand a function chosen to cancel the 
pole at q 2 = k 2 , leading to the result 

1 ^ fi m (k) - f lm (q)e ai<}2 - k2) m 



(2tt) 3 q 2 - k 2 

which is valid up to terms which are exponentially small in the volume. The exponential 
factors exp[a(q 2 — k 2 )} (with a > 0) are included so that the subtraction does not introduce 
ultraviolet divergences. Note that the factor of k l multiplying Y/ m is necessary so that 
the subtraction does not introduce a singularity at k — 0, and thus invalidate the use of 
eq. (J7J). 

From eq. (jlljl we immediately deduce the required summation formulae 4 



S lm (q) = 6 lfi ^ / 70 Y 00 + f lm {q)ZUq) , (12) 



3 Since all variables in this subsection correspond to the centre-of-mass frame, we do not introduce the 
symbol * which we use elsewhere in this paper to distinguish centre-of-mass variables from those in a 
moving frame. 

4 As above, this equality holds up to exponentially small corrections, but here, and in similar results 
below, this will not be stated explicitly. 
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where 

1 e a{q 2 ~k 2 ) . J3u e a(q 2 -k 2 ) 

Zi m {q) = T3 E -T^ kl 0) - 5 l>0 V j (27T)3 g2 _ k2 loo • (13) 

Here we have used the result that the integral on the right-hand side of eq. (jllj) vanishes 
for I as the integrand is rotationally invariant. This integral does not vanish, however, 
for I = 0, and must be included. Although the integrand in eq. (jllj) has no pole at k = q, in 
eq. (jT2j) we separate it into two terms each of which does have such a pole. For consistency, 
the two terms need to be regulated in the same way and the principal value prescription, 
denoted V, is a natural choice. 

The result in eq. (j!2j) was obtained for I — in ref. [6] [eq. (C.24)] using different methods. 
Our approach here allows a simple generalisation to all I, and, as will be shown in the next 
subsection, a generalisation to a moving frame. The form of (j!2j) is intuitively reasonable, 
since, when the summand is singular and thus rapidly varying between grid points, one 
expects the sum and integral to differ, with the difference proportional to coefficient of 
the singularity. In terms of the Poisson summation formula, eq. ©, the presence of the 
singularity requires one to keep the terms with / 7^ 0. 

Our result for the kinematical function Zi m (q), eq. (JT3J), appears to depend on the param- 
eter a. This, however, is misleading: the a dependence only enters in the exponentially 
small terms which we do not control. 5 This follows from the derivation itself; the key 
result, eq. (fTTj). is valid for any a, up to exponentially suppressed terms. It can also be 
seen directly by evaluating the dependence of Zi m on a, which has the advantage of show- 
ing the explicit form of the a dependence. For example, using the Poisson summation 
formula one can show that 



dZ 00 e c 



y 00 E ex pH-^ 2 /(4«)]. (14) 



da (47ra) 3 / 2 

The general expression for dZi m /da has the same exponentials multiplied by powers of 
L. This demonstrates that, if L is sent to infinity with q and a fixed in physical units, 
the result for Z\ m is independent of a. In practice, when working at fixed L, one should 
choose a sufficiently small that the a-dependent terms are negligible. The simplest choice 
is to send a — > + . 

This is one example of the more general result that Zi m should be independent of the ultra- 
violet regulator. We have checked this numerically by verifying the equality of our results 
(in the limit a — > + ) to those obtained in Refs. [3] using zeta-function regularization. 6 

The summation formulae presented above contain those given in ref. [2]. The latter are 
given in terms of the sum 



-E' 



L 3 - q 2 - k 2 



(15) 



5 Note that for / = the sum and integral in the definition of Z o both have an ultraviolet divergence 
proportional to a~ 3 / 2 , but these cancel leaving the exponentially suppressed dependence on a. 
6 The precise relation between our Z\ m and the Zi m of refs. [3,7] is given in eq. 
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with q 2 chosen to be one of the values of k 2 and the prime denotes that the terms with 
k 2 = q 2 are removed from the sum. The formulae in ref. [2] can readily be obtained by 
taking the limit \q\ —>■ \K\, where K is one of the values of k appearing in the sum in 
eq. (|8jl. This is shown explicitly for the case I = m = in appendix C of ref. [6]. 



2.2 Summation Formulae in a Moving Frame 

We now generalise the discussion to a moving frame, with total energy-momentum (E, P). 
As shown below in sec.|3J we require the summation formula for sums of the form 

k 

— * 

Note that, while the summation is over the moving frame momenta k = {2n /L)n } with n 
being a vector of integers, it is convenient to rewrite the summand in terms of the centre- 
of-mass momenta k* using the Lorentz transformation of eq. (jH7j) . It is also convenient, as 
will be apparent shortly, to pull the Jacobian wt/uj k out of the function f(k* ). 

We now proceed as in sec. 12.11 writing k* is terms of spherical polar coordinates, k* = 
(k* , 8* , <j)*) , and expanding f(k*) is spherical harmonics 

oo I 

/(£*) = E E fim(k*)k* l V^Y lm (e*, <!>*). (17) 

1=0 m=-l 

The factor of \/47r = l/Y 00 is introduced to simplify the subsequent expressions for I = 0. 
Again using eq. © we can write down the required summation formula 

77^ ~ k ^¥~ 2 k ^ Ylmie A } 



k 



d 3 k Ut flm(k*) ~ flm(q*)e a(q * 2 - k * 2) 



FV47TF, 



(2tt) 3 u k q *2_ k *2 

q *2_ k *2 k ™ Ylm ( 6 ' ^ ) • ( 19 ) 

The Jacobian factor oj%/co k corresponds to the change of integration variables from the 
laboratory-frame momenta k to the centre-of-mass frame momenta k*. The integral in 
eq. (j!9|) is the same as that appearing in the derivation of the centre-of-mass summation 
formulae (apart from the factor of \/47r), i.e. the same as the integral in eq. (jlljl . It is 
only non- vanishing for / = 0. 

The summation formulae above can be rewritten into a form convenient for subsequent 
manipulations. For I = we have 

1 Ljt e «(9* 2 -fc* 2 ) r d^k* e a(q* 2 -k* 2 ) 

C P (q* 2 )=cUq* 2 ) = — 9 *2_ fc *2 - V J J^3 q *2_ k *2 > ( 21 ) 
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while for I > the result is 




A- 



flm( k *) 

u>k q* 2 — k* 2 



(22) 



(23) 



A 



q* 2 — k* 2 



Note that we use c p as a shorthand for Cq Q . As in the centre-of-mass frame, the derivation 
given above shows that the quantities c[ m (q*) are independent of a up to exponentially 
small corrections in the volume. It is straightforward to use the expressions above to 
evaluate the cf m numerically. 

It is more difficult than in the centre-of-mass frame to use the Poisson summation formula 
to obtain simple analytic expressions for the a dependence of the cf m , since the terms with 
/ 7^ in eq. © involve a Fourier transform with respect to k and not k* . We have studied 
this issue numerically and have found that the dependence on a for even I is similar to 
that found in the centre-of-mass frame (see eq. (jl4j) ). For odd / on the other hand, the 
dependence on a appears to be linear, consistent with the presence of terms proportional 
to aexp(— mL) 7 . These differences however, do not affect the essential points necessary 
for control of finite volume effects. We repeat that for fixed a and physical parameters, 
the cf m are independent of a up to exponentially small corrections in the volume. In 
simulations performed at (large) fixed volume, a should be taken to be sufficiently small 
so as not to enhance numerically the terms which are formally exponentially small in the 
volume. This can be achieved for example, by taking the limit a — > 0+. 

In summary the summation formula for the sum in eq. ()16)) is 



3 Quantization Condition in Moving Frames 

The two-pion spectrum in finite volume can be determined from the exponential depen- 
dence of correlation functions of composite operators. We use 



(with time ordering implicit), where cr(x, t) is an interpolating operator for two-pion states 
and crp(t) is its spatial Fourier transform: 




(24) 



Cp{t) = (0\ap{t) <x f (0,0) |0) 



(25) 




(26) 



We will consider the two-pion correlation function in energy space 




(27) 



7 We have checked this in detail for P lying along an single lattice axis, and for I < 6. 
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In infinite volume Cp(E) has a two-pion cut in the s-plane (where s = E 2 — P 2 is the 
Mandlestam variable) starting at the branch-point at s = 4m 2 . The quantization of 
momenta in a finite volume means that the two-pion cut is replaced by a series of poles 
and the energies corresponding to these poles give the two-pion spectrum. Our aim in 
this section is to determine the position of these poles. We restrict the energy E to be in 
the range < s < 16m 2 , so that we can neglect finite volume corrections from four-pion 
or higher multiplicity intermediate states. The lower limit on s ensures that the energy 
is positive in all frames. 

Although we phrase our discussion in Minkowski space, we note that this same object 
may be obtained from the Euclidean space correlators calculated in lattice simulations 
by analytic continuation to imaginary Euclidean energy (which is the approach used in 
ref. [2]). 

The operator a(x, t) must be chosen to have overlap with whatever lattice two pion states 
we are interested in, but the details are not important here. We refer to ref. [7] for 
discussion of the symmetry group of two pion states for different choices of P. 

For simplicity, and to match the treatment in refs. [2,7], we will take the two pions to be 
indistinguishable. Although this restricts the scattering to even I, we will not display this 
constraint explicitly. This is because the final formulae can be shown straightforwardly 
to hold for all I in the case of distinguishable particles, although there are differences by 
factors of 2 at intermediate stages in the derivation. We will also assume that there is 
a symmetry, analogous to G-parity for physical pions, which forbids intermediate states 
with odd numbers of pions. We stress, however, that we make no assumptions about the 
form of the pion interaction, e.g. we do not rely on chiral perturbation theory or any 
other expansion scheme. Thus our results are completely general and apply to any two 
particle system. 

It will be helpful below to use the relations between the centre-of-mass variables and those 
in the moving frame. E and P are the total energy and momentum in the moving frame 
and E* = \JE 2 — P 2 = y/s is the corresponding total centre-of-mass energy. The relative 
velocity (J3) between the centre-of-mass and moving frames is given by 

P = P/E, (28) 

and E* and E are related by 

E = jE* where 7 = 1/^1 - (3 2 . (29) 

It is also useful to recall the definition of q*, the magnitude of the relative momentum in 
the centre-of-mass frame 

q* 2 = ^E* 2 -m 2 . (30) 

We stress that, in the following, the finite volume summations are always over the moving- 
frame momenta k = (2tt/L) n, with n a vector of integers. 

The correlation function Cp can be expressed in terms of the Bethe-Salpeter kernel K 
through the series shown in fig.d Since we choose E to lie below the four-pion threshold, 
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Figure 1: Diagrammatic expansion of the correlator Cp(E). Propagators are fully dressed 
and normalized to unity on shell. K is the amputated two-particle irreducible four-particle 
correlation function. The circles at the ends represent the operator a, renormalized by 
two factors of \fZ. 

there are no intermediate states with four or more pions and the finite-volume effects 
in K are exponentially suppressed [2,6]. The same is true of the dressed single particle 
propagators [1]. The only power-law volume corrections arise through the the two pion 
loops, and we now turn to an analysis of these corrections. 

3.1 Finite Volume Effects in the Generic Loop Summation 

The generic loop integration/summation appearing in fig. Qis of the form 

r = J_y- f dk o f(h,k) , , 

~ L 3 Y J 2tt (k 2 -m 2 + ie)((P-k) 2 -m 2 + ie) 1 } 

where k = (k , k) and P = (E, P) are four vectors and we have left out a factor of 
i 2 = — 1 from the propagators which will be accounted for later. The function / contains 
the energy-momentum dependence arising from the kernels on either "side" of the loop 
(with one of these kernels replaced by the matrix element of the operator o if the loop 
lies at one of the "ends" of bubble chains shown in fig. as well as that from the dressed 
propagators. Since E and P are fixed, / can be expressed as a function of the four- vector 
(ho, k) alone. Note that the dressing of the propagators does not effect the allowed values 
of the three momenta, which are constrained as usual by the finite volume. The only 
properties of / that we need are, first, that it has no singularities for real k (which holds 
given our kinematical constraint on E) , and, second, that its ultraviolet behaviour is such 
as to render the integral and sum convergent. 

The power-law volume corrections arise, as usual, because the integrand/summand is 
singular. To simplify the pole structure we first perform the ko integration. We choose to 
close the contour of integration so as to pick up the "particle" contribution from the first 
pole and the "anti-particle" contribution from the second: 

r ;J_yl /Kj) f(E + uj Pk ,k) \ 

I? Y \ ^{{E - u k ) 2 - u Pk ) + 2u; Pk ((E + u Pk ) 2 - ul) j ' 1 ] 

where 



oo k = y k 2 + m 2 and up k = y (P — k) 2 + m? . (33) 
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We have dropped the factors of is since what remains is a summation, and there is no 
obstruction to setting e = 0, as long as the poles do not coincide with the allowed values 
of k. This will be true in general, since the energies of interest are shifted by interactions 
away from those of two free particles. This argument does not hold in infinite volume, 
where we must perform an integral over k, and the factors of is must be kept. We return 
to this point below. 

For the kinematic region of interest, < E 2 — P 2 < 16m 2 , it is straightforward to show 
that the only singularity in / is the explicit pole in the first term inside the braces in 
eq. (|32|). which occurs at those values of E for which there is a term in the summation 
with uo k + ujpk = E. It is this singularity which leads to finite- volume corrections which 
decrease like powers of the volume. There are no other poles (either explicit or hidden 
in /) for E in this kinematical range. For the second term in eq. (|32|) we can therefore 
replace the summation by the corresponding integration and write 



where 



-i 



I = h + h, (34) 
1 y /Kj) and ( 35 ) 



d 3 k f(E + oo Pk ,k) 



(2vr) 3 2u Pk ((E + u Pk ) 2 -u 2 ) 



(36) 



To determine the finite-volume corrections we need to examine I\ in more detail. It is 
convenient to change variables in the summand to those in the centre-of-mass frame. In 
the moving frame k = (k\\,k±) is the momentum of an on-shell particle with mass m, 
where k\\ and k± are the components parallel and perpendicular to (3 respectively. The 
corresponding momentum and energy in the centre-of-mass frame, {kpk±) and u k , are 
given by 



fcjj = "f(k\\ — PuJk), k* ± = ki_ and uj k = \Jk* 2 + m? = ^{uJk — (3k\\) . (37) 
We now rewrite 1% in the form 

/ j 1 1 y 1 ,1 /*(£*) E* + 2u* k 

1 L 3 E* Y 2u k E* - 2u* k L 3 2E* ^ u k q* 2 - k* 2 Auj* k ' 1 ' 

where /* is the function / rewritten in terms of the centre-of-mass variables. Since cu k is 
a dependent variable, we choose to write /* as a function of k* alone. 

As expected, we see that the poles correspond to those terms for which io k = E*/2 or 
equivalently k* 2 = q* 2 . We have deliberately not rewritten the factor l/u> k in terms of 
the centre-of mass variables, so that the summand contains the same Jacobian as in the 
summation formulae (|2U|) and (|22j). Using these formulae, we find 
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where we recall that f* m is defined as in eq. (JT7J) and c^(g* 2 ) are defined in eqs. (}2Tj) and 

©• 

There is one further step we need to perform in order to write l\ as the infinite-volume 
result together with a correction, i.e. we should replace the principal-value integral in 
eq. ()39)1 by the corresponding integral with the Feynman is prescription in the propagator 
and a "delta-function" term: 

. 1 r d 3 k* /*(£*) E* + 2u* k 
1 ~~ ~ l 2E* J (2tt) 3 g *2 _ k* 2 + ie 4u? k 

+ -^EE /L(<f ) cL(g* 2 ) . (40) 

1=0 m=—l 

Note that the "delta-function" term picks out the I = part of /*. Observing that the 
first term in eq. (|4()|) is exactly the infinite volume expression for I\ in Minkowski space 
(after retracing the steps in the derivation above), we arrive at our final result, 

/ = Ioo + Ifv (41) 

^ (=0 m=—l ) 

where 1^ is the infinite volume result for the original loop integral, 1^ — Ii j00 + h- 

The quantity Ipy is the desired finite volume correction. Note that it involves the function 
/* with momentum argument k* = q*, so that both particles are on mass shell. Thus 
the infinite volume amplitudes contained in /* are evaluated with physical momenta, 
just as if one were picking out the cut in an infinite volume correlator. What is striking 
about this result is that q* is a momentum determined by the choice of E, and can take 
any positive value. In particular, q* is not constrained by the finiteness of the volume: 
it does not have to take on one of the (boosted) free particle values. This is a crucial 
result in the next section where it allows the finite volume corrections to be expressed 
in terms of physical scattering amplitudes, with the finite volume entering through the 
purely kinematic functions cf m . 



3.2 Quantization Condition 

In this section we use the result in eqs. (j41|) and (}4"2*|) to derive an expression for the two- 
pion energies in finite volume. We first derive the general, though formal, result, next 
work out the simplest case in which only the s-wave scattering amplitude is non-vanishing, 
and finally give the result assuming that the scattering amplitude vanishes for I > l max . 

Our approach is straightforward: we use eq. (J41j) for each of the loop summations in fig. ^ 
and then reorganize the series. If we keep in each loop then we recover the infinite 
volume correlator C'p(E). As already discussed, it does not contain the poles which we 
are seeking. The part of interest is the finite volume correction 

C F p v (E)=C p (E)-Cf(E). (43) 
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Figure 2: Contributions to the volume dependent part of the o correlator, Cp V . The 
notation is as in fig. ^except that the filled circles represent the full scattering amplitude, 
M, given by a geometric sum of any number of insertions of the kernel K, and the vertical 
dashed lines indicate that the on-shell finite volume part, Ifv-, nas been used for the loop 
integral. The quantities A and A' are defined in the text. 

This contains the desired poles, as well as cuts which must cancel those in Cp(E). Dia- 
grammatically, Cp V is obtained by keeping at least one insertion of Ifv- These contribu- 
tions are shown in fig.0 and lead to the following general result: 

Cp V (E) = -A' F A + A' F (zAf/2) F A+ ... (44) 

= -A' F A . (45) 

1 + iMF/2 K ' 

Here we have taken into account the factor of i 2 dropped from the loop in the previous 
section, as well as symmetry factors of 1/2 arising from the identical nature of the particles. 

The quantities in eqs. f)44j) and (|45jl act on the vector space of states of two pions with 
total energy E and momentum P. In particular, A and A' are given by 

A = out(7T7r ; E, P\a\0, 6)\0)Z n , A' = Z 7r (0|cr(0, 6)\mr ; E, P) in , (46) 

and involve a geometric series of contributions with any number of insertions of the kernel 
K, with the intermediate loops containing the infinite-volume expression 1^. The factor 
of ie in 1^ leads to the presence of in- and out-states as shown. M is the on-shell scattering 
amplitude, obtained by summing a similar geometric series. And, finally, F is a kinematic 
factor determined by the form of Ifv-i which, in particular, accounts for the fact that the 
"cut" can have non-vanishing angular momentum. 

To give these quantities a concrete representation we change to centre-of-mass variables, 
since we can then use the partial wave basis. In particular, M is diagonal in this basis, 
with matrix elements 

M —A A 16tt£* (exp[2t5 h (q*)] - 1) 

q* 2l 

where 5i is the phase shift in the Z-th partial wave. For Ipv we have already changed to 
centre-of-mass variables in the result (J42|) . from which it follows, using the definition of 
fim [ see ec f (E3)L an d t ne completeness of the spherical harmonics, that 



h,mi;l2,m2 



q 

W 7 f\ 7 V •/ 
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where the implicit arguments of the spherical harmonics are (9*, 0*), and the asterixes on 
the spherical harmonics indicate complex conjugation. The integral can be expressed in 
terms of Wigner 3-j symbols, but an explicit expression will not be needed. 

As for the unknown matrix elements A and A', the change of variables leads to no sim- 
plification, but no complication either. They are simply unknown vectors A\ m and A' lm . 
Note that, unlike M, the matrix elements A and A' need not be Lorentz invariant, since 
the operator a can involve spatially separated pion fields. But this lack of invariance is 
of no consequence for our argument, since the role of A and A' is simply to provide the 
coupling of the external operators to the on-shell two pion states. In other words, we are 
interested here in the positions of the poles, and not their residues. 

The positions of the poles in Cp V are determined by the factor which lies between A' and 
A, namely F(l + iMF/2)~ 1 . F itself inherits from cf m poles at all energies corresponding 
to two free pions both having allowed momenta in the moving frame, as is clear from 
eqs. (12 1| 123 (1 . Turning on the interaction will, however, shift these poles, and indeed one 
can see that they cancel between numerator and denominator if M ^ 0. The poles are 
shifted to the positions determined by requirement that 1 + iMF/2 have zero eigenvalue, 
or equivalently by the formal quantization condition: 

det(l + iMF/2) = 0. (50) 

This is the main result of this section. It demonstrates that, for any P, the finite volume 
energy shifts depend on interactions only through the infinite volume scattering amplitude, 
as was shown by Luscher in ref. [3] for P = 0. The result ()50|) also gives a clear separation 
between the dynamical effects (the scattering amplitude M, independent of the volume) 
and the kinematical effects of finite volume (the matrix F, dependent on the volume, but 
independent of the scattering amplitude). 

In the derivation of ref. [3] it was important that the range of the effective non-relativistic 
potential, R, satisfy L/2 > R, aside from an exponentially falling tail. Without this 
condition there is no "outside" region where the wavefunction is composed of free waves 
with phase determined by the phase shift. Thus for L/2 < R one would not expect the 
energy shift to depend only on the phase shifts, i.e. on M alone. By contrast, our result, 
eq. ([50)1 . appears to imply that the energy shifts depend only on M without any proviso. 
This is, however, misleading. Our entire approach relies on the assumption that sums of 
non-singular functions in momentum space can be well approximated by integrals. When 
the range of the potential is larger than the box size, the spread of its Fourier transform 
(Ak ~ l/L) is smaller than the spacing between the discrete finite volume momenta. 
Consequently, the starting assumption of our derivation fails, and there is no conflict 
with the result of ref. [3]. 

To make the result (|5Uj) into a practical tool we need, following refs. [3,6,7], to truncate 
the number of partial waves that contribute. In other words, we assume that 5i = for 
I > ^max- In the following subsection we discuss the simplest case, Z max = 0, and then 
return to the general case where < / max < oo. 
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3.2.1 Quantization Condition: s-wave Interactions 



We assume that A and A' couple the external operators to two pions with an s-wave 
component, and that only the s-wave part of the scattering amplitude is no n- vanishing. 
Note that although F couples the s-wave to other waves [through the cf m terms with 
/ 7^ in F FV , eq. (J32J)], these couplings would have to contribute through terms of the 
form F M F, F M F M F, etc. in eq. (144)1 but these vanish since M = for the higher 
waves. Thus the geometric series leading to the pole occurs only in the s-wave, and can 
be treated in isolation. We discuss this point more thoroughly in the general context of 
the next subsection. 

Since the problem is now one-dimensional, the quantization condition is simply 



where 



Ms 



1 + iM s F s /2 = , 
_ lQnE* (e™(^ - 1 



(51) 
(52) 

8«E* 2E*~ - > ™ 

We note that the quantization condition ()51j) has been obtained in ref. [12] using a similar 
method to ours in the context of a non-relativistic effective theory of nucleon-nucleon 
scattering. After rearrangement, the quantization condition can be written in the form 



q* 2% 

with 5 the s-wave phase shift, and F s is the / = component of F: 



00:00 



tan [<5(<?*)] = — tan 4> p (q 



where tan 



q*_ 

4:71 



(54) 



Thus the finite- volume two-pion states on a cubic lattice of size L 3 are those with energies 
(and hence values of q*) which satisfy eq. (j54j) . For the particular case of P = one 
immediately recovers the Luscher condition [3] . 

We note in passing that for external operators in the moving frame having no centre-of- 
mass s-wave component, the allowed finite volume energies are those of two free particles. 



3.2.2 Quantization Condition: general interactions 

In this subsection we assume that the phase shifts, and consequently scattering ampli- 
tudes, are small enough that they can be ignored for I > / max . In other words, we assume 
that 

^l,m;l,m = if I > l max ■ (55) 

As we now explain, this turns out to imply that one can truncate the entire problem 
to the subspace with I < Z max . In particular the determinant of the full matrix (1 + 
iMF/2)i ltmi -i 2tm2 vanishes if the determinant of the block with < ^max vanishes. 
Thus the quantization condition becomes 



det(l + iMF/2) 



= 0. (56) 

h,h<i max 
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The condition (foTj) of the previous section is the simplest example (/ max = 0) of this 
general result. 

The subtlety in establishing this result arises from the fact that F has non-vanishing 
elements connecting li < / max and Z 2 > Lx an d vice-versa. To see that these "block-off- 
diagonal" matrix elements are irrelevant, it is useful to introduce a projection operator, 
P, onto the subspace with I < / max . Then the assumed form of the scattering amplitude, 
eq. (|55|) . can be written as M = PMP. Using this property, and denoting the truncation 
of F to the / < / max block by Fp = PFP, one can manipulate the result for the finite- 
volume part of the correlation function into the form 



cy(E) = -A' 



(i/2)FM 



P + iMFp/2' 



A. (57) 



This shows that poles occur when the truncated matrix P + iMFp/2 has zero eigenvalues, 
which is equivalent to the condition (JoTlj) . 

A consequence of this argument is that we only need to determine F in the truncated 
basis. Using the explicit form in eqs.pS|) and (J2Hj), we see that if h.h < ^max, then the 
angular integral in F FV vanishes for I > Z1 + /2. Thus only a finite number of the quantities 
cf m are needed to calculate the truncated determinant. 

In order to compare our result to that of refs. [3, 7], it is convenient to pull out the 
volume- independent first term from F in eq. (J48j) . Then one can rewrite the quantization 
condition as 



det(cos£ — sin 5 F 



= 0, (5* 

h,h<l max 



where cos 5 is a diagonal matrix with entries cos 5^ on the diagonal, and similarly for 
sin S. This condition has the same form as eq. (100) of ref. [7], and agrees if our matrix 
F FV equals that denoted M. there [and by Ai RG in eq. ()59|) below]. We have checked, 
using the relation between our cf m and the zeta-functions Zf m of ref. [7] given in eq. (J73)l 
below, that in fact, 

The extra factors of i cancel in the determinant in eq. ([58)1 . because cos 5 and sin 5 are 
diagonal matrices. Thus our quantization condition agrees with that of ref. [7] . 



4 Comparison of Quantization Condition with that 
of Rummukainen and Gottlieb 

In this section we show that the kinematical functions cf m (q* 2 ), defined in eqs. (|21j) and 
(J23|) . can be rewritten in the form given in ref. [7]. As already noted in the previous section, 
this implies the equivalence of the quantization conditions. To compare the results requires 
a small kinematical exercise, i.e. the rewriting of the factor {uj^/ujk) / '{q* — k* ), which 
appears in the definition of cf m (q* 2 ) in terms of the momentum r of eq. ©• We need this 
in the vicinity of the pole, i.e. for q* 2 w k* 2 . We first recall that, for a given choice of k, 
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setting q* 2 = k* 2 corresponds to a kinematical situation with both pions on-shell. This 
follows from the derivation in the previous section, in which the pole in q* arises when 
both intermediate propagators are on shell. The particular consequence we need here is 
that u k = E*/2 at the pole, since the two on-shell particles share the energy equally in 
the CM frame. This result can be seen directly from 

u* k 2 - (E*/2) 2 = (k* 2 + m 2 ) -(q* 2 + m 2 ) = k* 2 - q* 2 . (60) 

An important corollary is that f = k* at the pole. To see this we recall from ref. [7] that 

r y =7- 1 (/c|| -P/2), r ± = k ± = k* ± . (61) 

Using the inverse Lorentz transform, k\\ = j{k*\ + (3u k ), together with E* = E/j and 
(3 = P/E, one finds 

ry = kl + (3u* k - P/(2 T ) = fcj + P(u* k - E*/2) . (62) 

Since u k = E*/2 at the pole, this establishes the desired result. Thus setting q* 2 = k* 2 
implies also q* 2 = r 2 , and it follows that the positions of the poles in q* in the sum 
contained in c p , eq. (f2*Tjl. are the same as those in the sum of ref. [7], eq. (J1J). This is 
clearly a necessary condition for the equivalence of the results. 

It is now straightforward to determine the residue of the pole when written in terms of 
r 2 rather than k* 2 . Working to linear order in the difference E* /2 — u k we have 

/?* 2 I A. — i i* 2 n* 2 — h* 2 

Er/2-utM - /4 Uk = q - — . (63) 

2w* 2u* k v ; 

Using this result in eq. ((02} yields 

k* 2 - r 2 » 2(3k\{E*/2 - u* k ) « (3k\{q* 2 - k* 2 )/u* k . (64) 

Thus we find 

q* 2 -r 2 = (q* 2 -k* 2 )(l+(3kl/u* k )+0[(q* 2 -k* 2 ) 2 ] (65) 

= (g* 2 - k* 2 )(u k / 7 u* k ) + O [Of 2 - F 2 ) 2 ] , (66) 

where in the second step we have used U/- = 7(a>£ + fiku). 

Using this result we can rewrite our kinematical functions in terms of r 2 . We begin with 

c p (q* 2 ): 

p 1 u k e Q (9* 2 ~ fc * 2 ) r d^k u* (f(i* 2 ~ k * 2 ) 

° ( Q ^ = L3 E — q *2_ k *2 ~ V J J^3 — q *2_ k *2 ( 6T ) 

£ \-^-Z 3-+»(*) "Mt^s £- J _- 5 . + ,(*) ■ (68) 



7L 3 4^ [ g* 2 — r 2 J 7 7 (27r) 3 [ g* 2 — r 



Here g is a non-singular function of k, which arises from the 0[(q* 2 — A;* 2 ) 2 ] terms in 
eq. ()66)1 . with contributions both from the pole and the exponential. It has an implicit 
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dependence on q*. Its form is not required, however, since the sum and integral over g(k) 
cancel (up to exponentially small terms in the volume) as a consequence of the Poisson 
summation formula. Thus we find that 



e a(q* 2 ~r 2 ) I r d 3 k e^'^ 



c P («* 2 ) = -PT—I ' (69) 



7L 3 ~f q* 2 — r 2 7 J (27r) 3 q* 2 — r 

With this expression for c p (g* 2 ), the quantization condition in eq. (|54)l is manifestly equiv- 
alent to that of ref. [7]. 8 

The equivalence can be extended straightforwardly to higher partial waves. The key point 
is that the vectors k* and r are equal at the pole, as established above. Thus one can 
replace angular dependence on k* with the same dependence on r, with the difference 
being non-singular. The argument goes as follows (assuming / > 0): 

- is? ZiT^^ 1 "™ 

A: 

1 f e a (<?* 2 - r2 ) _ 1 

Eh] T r* l VbrY lm {6r,<f>r)+gim{k)\ 



1 r d 3 k ( f a (i* 2 - r2 ) , , ^1 

~j V Jj^\-^^ r ^ Y ^^+9Uk)\ (72) 

1 e »(q* 2 -r 2 ) 

= -F^T,— T r* l V4^Y lm (6 r ,<f> r ). (73) 

In the first step we have added an integral that vanishes, but is nevertheless useful for 
subsequent manipulations. In the next step we have rewritten the expression in terms 
of r [using spherical components (r,9 r ,<f> r )], leading to the presence of the non-singular 
function gi m (k). To make the final step, we note that the sum and integral over gi m 
cancel, and that the remaining integral over k, which can be rewritten as an integral over 
r , vanishes by rotational symmetry. The result is an expression given in ref. [7] , although 
regularized differently. The precise relation is 



cL(q* 2 ) = f^V Vjl; Of V2vr) 2 ] , (74) 



7 L 3 \L 

where as usual the equality holds up to terms which vanish exponentially with the box 
size. As already noted in sec. 13.2.21 with this relation, our general quantization condition 
eq. (J5fi)l agrees with that obtained in ref. [7]. 

8 We use a different regularization, but the result should be independent of the ultraviolet regulator, and 
we have checked that our result agrees numerically with that obtained using the analytic regularization 
of ref. [7], up to terms which vanish exponentially as L — > oo. 
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5 Quantization Condition in Perturbation Theory 



In this section we test the validity of the quantization condition in eq. ()54|) in perturbation 
theory. We assume an interaction potential A7r 4 /4! and work to linear order in A. Specif- 
ically, we take the expressions for the s-wave phase-shift and two pion energy evaluated 
at 0(A) and verify that they satisfy eq. ()54|). 

The 7T7T phase-shift at 0(A) is readily calculable and is given by 

4 C> = -isf- (75) 

We stress that, for this exercise, we are considering indistiguishable "pions" without an 
isospin degree of freedom. 

We now sketch the derivation of the two-pion energy at one-loop order and then in sec. 15. 21 
check that they satisfy the quantization condition. 



5.1 Evaluation of the Two-Pion Energy in Perturbation Theory 

For the perturbative evaluation of the energy eigenvalues in a finite volume, we follow 
the steps in the derivation presented in section 4 of ref. [8], generalised to an arbitrary 
frame. We refer the reader to ref. [8] for a more complete description of the calculation. 
The energy eigenvalues are obtained by evaluating the time dependence of the Euclidean 
correlation function: 9 

= (0|7r^(f)7r^(t)7r(0) 2 |0), (76) 

where 

7Tf(t)= / d 3 xn(x,t)e iq ^ (77) 

is the Fourier transform of the pion field, which creates a pion with momentum q and 
destroys a pion with momentum —q. The momenta qi must take one of the values allowed 
by the box size, but are otherwise completely general. In particular, P = q\ + q 2 need not 
vanish. 

The two-pion energy can be deduced from the exponential dependence of the correlation 
function on t. Without loss of generality, we take t > 0. At tree level C^ 2 (t) can be 
represented by the diagram in fig.|3](a), and one obtains 



p—Ult p -u 2 t 

Cf\ (t) = 2- , (7c 

9192 v ' 2cji 2uo 



where ujf = q^ + m 2 , and the factor of 2 arises since there are two possible contractions. 
The two-pion energy is thus ui\ + io%. 



9 In ref. [8] the discussion considered the K — > tttt amplitude, but for the calculation of two pion 
energies we do not need the initial kaon. 
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(a) (b) 



Figure 3: Contributions to the correlation function defined in eq. (f7B|): (a) tree-level 
diagram; (b) one- loop diagram giving the energy shift at 0(A). 

At one-loop order, only the diagram in fig.[3j(b) gives rise to a shift in the two-pion energy. 
Following ref. [8], the contribution linear in t is seen to be 

t ^ - U 't> -uj' 2 t' - Wl (t-f) - U2 (t-f) 



Si'. 

on — shell 



Here g^' and q 2 are the loop momenta (with the corresponding energies), and the 
sum over q/ is over all "on-shell" intermediate states. These are the choices satisfying 
Qi + Qz = P ari( i uj[ + uj' 2 = uji + u>2- The number of such choices is denoted v. In the CM 
frame, this constraint implies u>[ = u 2 = uj\ = u 2 , so that all contributions are related 
by rotations, and all come with the same kinematical factor. For F ^ 0, however, the 
on-shell constraint can be satisfied in two ways: (i) those related by a lab frame symmetry 
transformation, for which u[ = u)\ or u 2 = lo 2 \ and (ii) accidental choices of q\' not related 
by a lab frame symmetry transformation, for which uj[ equals neither uj\ or u 2 , although 
u[ + uj 2 = uj\ + u 2 . The latter are absent in general, since they require special values of the 
pion masses. 10 They are also absent for the simplest choices of momenta, e.g. q\ = (1, 0, 0) 
and q 2 = (0, 0, 0) (in units of 2tt/L), for which v = 2 for all pion masses. 

In the following, we assume that there are no accidental on-shell configurations. In this 
case, the kinematical factor is common to all terms in the sum over t' and one finds that 
the term linear in t in is 

Combining this with the result for one finds, the energy shift AE = E — uji — uj 2 to 
be 

^ = jh^- (si) 

L A 8a>iu;2 

Note that this holds for any allowed choice of uj{. 

5.2 Verification of the Quantization Condition 

We recall that the quantization condition is 47r tan(5)/g* = l/c p (q* 2 ), since we only 
have s-wave scattering. Verifying this turns out to be simpler if one uses the form of 

10 One example is qi = (2,4,0), q 2 = (-2,1,0), = (1,5,0) and q 2 ' = (-1,0,0) with m 2 = 209/16 
(all quantities in units of 2n/L). 
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c p (q* 2 ) given in eq. ([69)1 . i.e. the form obtained in ref. [7]. Schematically this form is 
c p = 7 _1 L~ 3 J^kiq* 2 ~ r 2 )^ 1 . At leading order this condition is satisfied almost trivially. 
The allowed energies are those of two free pions each with momentum consistent with the 
boundary conditions. For any such choice we have q* 2 = k* 2 = r 2 , as discussed in sec. EJ 
Thus cp is infinite, requiring tan S = 0, which is consistent with the leading order result 
that 5 = 0. 

At next order we consider energies lying close to one of free pion values: E = U1 + UO2 + AE 



where LO{ = y q 2 + m 2 , with ^ being the allowed finite- volume momenta. The energy shift 
AE is of 0(A). Thus in each of the v terms in the summation which at leading order 
had q* 2 = r 2 the denominators q* 2 — r 2 are now of 0(A), and hence Cp is of 0(1/A). 
All other terms in the summation, as well as the principal part integral required for 
regularization, are of 0(1) and thus subleading. Keeping only the leading terms, the 
quantization condition becomes 



q* 7L 3 d(q* 2 -r 2 ^ 



AE 



Aw v dE 



(82) 



Here we have made use of the result, to be established below, that all v relevant terms 
give the same contribution to c p . 

The first part of the derivative is readily evaluated. Recalling that 4(g* 2 + m 2 ) = E* 2 = 
E 2 — P 2 , and noting that the derivative is to be performed at fixed P and m, we have 



dq* 2 _ E 
~~dE ~ 2~ ' 

Next, using the definition (jHTjl for r, and noting the k is held fixed, we have 



33) 



dr 2 drl d-v' 2 IP 2 



dE dE y 11 1 ' dE v 11 E 3 
Combining (JH3|) and we fi n d 

d(q* 2 -r 2 ) _ E A - (2fcjj - PfP 2 
dE ~ 2EP 



55) 



This expression simplifies when evaluated for two on-shell pions (i.e. for q* 2 = r 2 ). The 
four-vectors of the two pions in the lab frame are then 

<g = (u u k), q% = (u 2 , P - k) . (86) 

Observing that 

= (qi + g 2 r (gi - g 2 ) M = E(ui - uj 2 ) - (2k\\ - P)P , (87) 
the derivative at the pole is readily seen to be 



d(q* 2 — r 
dE 



E 2 — [uji — U2) 2 2ujiuj 2 
2E ~ E 
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Note that this result is the same for all v contributions to the sum in c p . 
Inserting this result in the quantization condition and rearranging gives 

n* A P J 3 

4 = TSF — 8 ^ 2 ' (89) 

This is consistent with the results for S and AE given respectively in eqs. ()75|) and (JSTj) 
above, confirming the validity of the quantization condition at 0(A) in perturbation the- 
ory. 



6 Finite- Volume Effects in Matrix Elements 

In this section we generalise the Lellouch-Luscher formula [5] for the finite-volume effects 
in K — > 7T7T matrix elements from the centre-of-mass frame to a moving frame in which 
the two-pions have total momentum P. We find the following relation between the matrix 
elements of the effective Weak Hamiltonian density, Tiw, in infinite and finite volumes: 

\A W \ 2 = 8nV 2 ^ {<%*) + \Mw\ 2 , (90) 

where the ' represents the derivative with respect to q*, 5 is the s-wave phase shift, (f) p the 
kinematic function defined in eq. (|54|) . V = L 3 is the spatial volume, and the K — ► titt 
matrix elements in infinite and finite volumes are respectively 

A w = ooin^E^lHw^lK;?)^ (91) 

and 

M w = v (Tnr,E,P\Hw{0) \K\P) V (92) 

(as indicated by the subscripts). The infinite volume states have the usual relativistic 
normalizations, while the finite volume states are normalized to unity. Note that for this 
formula we are assuming that 7-^(0) does not insert energy or momentum. In particular, 
for both infinite and finite volume the energy of the two-pion state equals that of the 
initial kaon, E 2 = m 2 K + P 2 . We comment on the generalisation below. 

The result (J9(J|) is valid, as usual, up to corrections vanishing exponentially with L. It 
also assumes that one can neglect the scattering amplitude in the I = 2,4,... waves at 
E* = rriK- Note that this is a stronger assumption for P^O than for P = 0. In the latter 
case, the first higher partial wave which contributes has / = 4, due to the cubic symmetry. 
For a moving frame, however, the symmetry is reduced and the leading contribution comes 
from I = 2. Technically this arises because oc ^ 0, as shown explicitly in ref. [7]. 

In the centre-of-mass frame the original derivation of the formula in eq. ()90|) proceeded 
by using degenerate perturbation theory [5], in which the volume was chosen so that, in 
the absence of weak interactions, the kaon and the two-pion states are degenerate. This 
degeneracy is then broken by the weak interactions. In the derivation of ref. [6] the factor 
relating the infinite and finite volume matrix elements was interpreted in terms of the 
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Figure 4: 0(A) contribution to the tttt scattering amplitude (see text). The thick line rep- 
resents the propagator of a kaon, which is off-shell by an amount of 0(A). A w represents 
the infinite- volume K — ► tttt matrix element. 



density of states of the two-pion system. We now derive eq. (|90j) using both approaches 
in turn. For simplicity, we continue to treat the pions as identical, ignoring the isospin 
degree of freedom. The final result applies, however, separately for each isospin state 
(using the corresponding phase shift). 



6.1 Derivation using degenerate perturbation theory 

In this section we generalise the derivation of ref. [5] for the finite- volume effects in matrix 
elements to a moving frame. We start by considering the case with no weak interactions 
(so that K — > tttt decays are not possible) and take the volume such that the energy of 
a kaon with momentum P is equal to the energy of one of the finite-volume two-pion 
states with the same total momentum. 11 The weak interactions are now introduced as 
an (arbitrarily weak) perturbation which breaks the degeneracy between the kaon and 
two-pion states. The energy eigenvalues are then 

E = Eq±\V\M w \, (93) 

where Eq = P 2 + m 2 K and A is a parameter introduced to facilitate the counting of orders 
of perturbation theory and is set to 1 in the final results. The factor of V is present 
because in eq. (I92|) we have defined My/ to be the matrix element of the Hamiltonian 
density. 

We can also calculate the shifts in (|9*H|) by taking into account the effect of the weak 
interactions on the phase shift, which then impacts the finite- volume energies through 
the quantization condition. The only contribution to the scattering amplitude (and thus 
phase shift) linear in A is that from the diagram shown in fig.Hl [5]. The energy and 
momentum of the external tttt states are chosen to be (E, P), where E is one of the finite 
volume energies in eq. ()93|). The kaon propagator is thus off shell by 0(A), and when 
combined with the two Kirn vertices, each of 0(A), the result is a contribution to the 
s-wave scattering amplitude of first order in perturbation theory: 



1-4 



|2 



^•^mm^'- < 94 > 



11 We are considering here only the two-pion states with an s-wave component, so that their energy is 
determined by the quantization condition eq. Ij54(l using the s-wave scattering amplitude. These states go 
over to pure s-wave states as the volume is sent to infinity. 
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where 5 represents the unperturbed phase-shift. Using eq. (|52j) we obtain the expression 
for 6, the perturbed phase-shift (i.e. the phase-shift including 0(A) terms): 

w = w*m m ?mh\ - (95) 

In order to obtain eq. (fHU|) we 

1. use the quantization condition to replace S(q*) by nn — (p p (q*)) 



2. expand <fi p {q*) and S(q*) around the unperturbed point q^ where Eq 2 — m 2 K = 
4go 2 + m 2 , 

<p p (q*) ~ <p p (q*) + Aq* <f> p '(q*) and 5(q*) ~ % *) + Ag* 5'(g Q *) (96) 

with 

An* = 4 

4g* 

3. Note that p (gg) = n7r — ^(<7o) anc ^ ec L ua t e the 0(A) terms on both sides of eq. 



XE 

Aq* = ±— V\M W \; (97) 



6.2 Derivation of eq. ( 1901) using the density of states 

In this section we derive eq. (190 jl by generalising the arguments of ref. [6]. Let a(x,t) be 
an interpolating operator which can create two-pion states from the vacuum and consider 
the two point correlation function in finite volume: 



C%(t) = I d 3 xe iP ^(0|a(x,t)a(0,0) \0) v (98) 
= V J2 I ( I °"(°) \^;P,n}\ 2 e~ Ent (99) 

n 

— > V ( dE p v (E) | ( | ct(0, 0) | tttt; E, P) v \ 2 e~ Et . (100) 

In eq. (IK)Ojl pv{E) is the density of states and can be evaluated using the quantization 
condition in eq. (|54jl . This condition can be rewritten as 

5(q*) + <p P (ql=nn, (101) 

with n an integer. We argue below that, for large enough L, 5+<f) p increases monotonically, 
so that n counts the number of states. Given this, the density of states is given by 

^ = % = ^M"")+rW)}. (102) 

The monotonicity of 5 + <p p follows because, as explained below, <p p is a monotonically 
increasing function of q* whose derivative is proportional to L. By contrast the dependence 
of 8 on q* is volume independent. Thus, although 5 need not be an increasing function of 
q*, its dependence will be overwhelmed by that of P for large enough L. 
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The behaviour of <p p can be seen most easily using the form obtained in ref . [7] : 



<j> p (q*) = tan- 1 [q*/(Anc p )} , c p = (^L 3 )- 1 ]>>* 2 - r 2 )' 1 . (103) 

k 

Aside from the factor of 7 _1 , c p decreases monotonically between the poles which occur 
whenever there is a free two-pion state in volume L 3 with total momentum P. As L 
increases, the spacing in q* between such states decreases, so the derivative c P/ becomes 
more negative, and overwhelms any positive contribution to the derivative from 7 _1 . This 
behaviour of c p implies that <p p increases monotonically at large enough q*. Furthermore, 
since (fi p increases by n for each transit between poles, whose spacing ~ 1/L, the derivative 
(p Pl is proportional to L. 

We now return to the main line of the argument. The infinite- volume correlation function 
corresponding to is given by 

C ™ {t) = vt^J dE -^\(°\°&°)\^ P > E )oo\ 2 e- Et , (104) 

where E* 2 = Aq* 2 + m 2 = E 2 — P 2 as above. Comparing eqs. (jlOOj) and ()104j) we obtain 
the correspondence between the infinite-volume and finite-volume kets: 



I nn; P, E ) co <£> 4n\ VE * Pv ^ | tctc; P,E) v . (105) 
V q* 

Similarly, by considering correlation functions of interpolating operators for single particle 
states, we find: 

\K;P) OQ ^V2EV\K;P) v , (106) 

with E 2 = P 2 + m 2 K . Eq. (|106|) is simply the relation between the wavefunctions of states 
normalized relativistically (2E particles per unit volume) and non-relativistically in a 
finite- volume (1 particle in volume V). 

Combining eqs. ()105|) and (jl06|) we obtain the required relation between the matrix ele- 
ments in infinite and finite- volume in eq. ()90|) . 

In eq. (|90j) the matrix elements correspond to physical processes and energy and momen- 
tum are conserved. However, eqs. (jl05|) and (|106J) can be applied more generally. One 
such application is to the chiral extrapolation of the matrix elements computed in lattice 
simulations. Since these are currently performed at unphysically large masses for the u 
and d quarks, the results have to be extrapolated to the physical point. Chiral pertur- 
bation theory (xPT) can be useful in guiding this extrapolation. For example, in ref. [8] 
it was proposed to fix the low-energy constants which appear in x?T at next-to-leading 
order by computing K — > tttt matrix elements with the kaon and one of the final state 
pions at rest and varying the momentum of the second pion. Equations (|105j) and (|106j) . 
with different energies and momenta for the kaon and two-pion system, can be combined 
to determine the finite-volume effects also in such cases. 
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7 Summary and Conclusions 



In this paper, we have provided a field theoretic derivation of the finite volume energy shift 
for two hadron states in a moving frame, confirming the result obtained by ref. [7] using 
a relativistic quantum mechanical approach. We have also determined the finite-volume 
corrections in matrix elements of local composite operators with an initial and/or final 
state consisting of two hadrons, thus generalising the Lellouch-Luscher factor to moving 
frames. The path is therefore now open to numerical studies of two pion energies and 
matrix elements in a moving frame. As described in the introduction, there are a number 
of important applications which depend upon, or are greatly simplified by, working with 
P ^ as well as P = 0. 

Our work provides three new results. The first is technical: we give a simple derivation of 
the summation formulae needed to determine finite- volume effects. It has the advantages 
of conceptual simplicity, and that it can be readily applied in any frame. In addition, it 
provides a straightforward way of calculating the required regularized sums numerically. 
It is of comparable numerical efficiency to the method used in refs. [3] and [7], based on 
analytic regularization. 

Our second new result is a method to determine the energies of finite volume states, to 
all orders in perturbation theory in 1/L, based entirely in field theory. This is in contrast 
to the method developed in refs. [2, 3] for the centre-of-mass frame. There it is shown 
that, in the centre-of-mass frame, the two-particle energies in a general relativistic field 
theory are related to those of an auxiliary non-relativistic quantum mechanical (NRQM) 
system in the same finite volume, with an energy-dependent potential chosen to reproduce 
the relativistic scattering phase shifts. The analysis is then carried out in the NRQM 
system. Perturbation theory in 1/L is developed in ref. [2], with ref. [3] providing the 
generalisation to all orders in 1/L. The use of a NRQM theory is a technical device, 
and does not require the underlying two-particle system to be non-relativistic. It reflects 
the fact that the dominant finite volume effects come from nearly on-shell two-particle 
intermediate states. Nevertheless, we believe that it is conceptually simpler to avoid the 
need for an auxiliary NRQM theory, and this is what our approach allows. In addition, 
we obtain the general result with multiple non-vanishing partial wave amplitudes in a 
very straightforward fashion. 

The particular advantage of our approach, however, is that it is readily generalisable to a 
moving frame. The general formalism remains unchanged, but the required kinematical 
functions, which arise from the relation between two-particle loop momentum sums and 
integrals, depend on P. The generalisation of the method using the auxiliary NRQM 
theory is less obvious, in part because of the mismatch between relativistic and non- 
relativistic energies for P ^ 0. This is finessed in ref. [7] by the use of a relativistic two 
particle wave function, satisfying the Klein-Gordon equation. The wavefunction is subject 
to periodic boundary conditions in the moving frame, and the key step is to convert these 
to boundary conditions in centre-of-mass coordinates. One would expect the latter to 
involve unequal times because of the Lorentz boost, but the authors of ref. [7] argue 
that this is not the case. Their argument is based on the observation that outside of 
the interaction region the wave function is that of two free particles and can readily be 
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shown to be independent of the relative time in the centre-of-mass frame. From this they 
conclude that the boundary conditions can be applied at equal times also in the centre-of- 
mass frame, with the effect of the boost causing an elongation of the box in the direction 
of P. We were left with the question as to whether this procedure correctly corresponds to 
periodic boundary conditions in the moving frame (including inside the interaction region) 
and this question partially motivated us to develop the alternative approach presented 
here. Since our final results for the quantization condition can be brought into the same 
form as those of ref. [7], we have indirectly confirmed their assumptions. 

As a check of our result (and thus also that of ref. [7]) we have calculated the two-pion 
energies to leading non-trivial order in the perturbative expansion in a A0 4 theory, and 
confirmed the quantization formula in a general moving frame. 

Our final new result is the generalization to moving frames of the Lellouch-Luscher factor 
which contains the finite-volume effects in K — > tttt matrix elements. This turns out to 
require only a simple kinematical change: the non-trivial effects of working in a mov- 
ing frame are already incorporated in the functions c[ m that determine the two-particle 
energies in finite volumes. 
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